Microscopic origin of local moments in a zinc-doped high-T c superconductor 



o 
o 

O 

Q 

o 

m 



X. L. Qi and Z. Y. Weng 
Center for Advanced Study, Tsinghua University, Beijing 100084, China 
(Dated: February 2, 2008) 

The formation of a local moment around a zinc impurity in the high-T c cuprate superconductors is 
studied within the framework of the bosonic resonating- valence-bond (RVB) description of the t — J 
model. A topological origin of the local moment has been shown based on the phase string effect 
in the bosonic RVB theory. It is found that such an S = 1/2 moment distributes near the zinc in a 
form of staggered magnetic moments at the copper sites. The corresponding magnetic properties, 
including NMR spin relaxation rate, uniform spin susceptibility, and dynamic spin susceptibility, 
etc., calculated based on the theory, are consistent with the experimental measurements. Our work 
suggests that the zinc substitution in the cuprates provide an important experimental evidence for 
the RVB nature of local physics in the original (zinc free) state. 

PACS numbers: 74.20.Mn,74.25.Ha,75.20.Hr 
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I. INTRODUCTION 

The Zn substitution of the in-plane Cu 2+ ions in the 
high-r c cuprates introduces some novel properties to the 
system. Although a Zn 2+ ion can be considered as a 
nonmagnetic impurity in the Cu02 plane, strong mag- 
netic signatures have been detected at the surrounding 
Cu sites. The NMR and NQR experimento 1 ' 2 ' 3 ' 4 ' 5 have 
shown the presence of staggered antiferromagnetic (AF) 
moments near the Zn site, whose sum behaves like an 
5=1/2 magnetic moment as indicated by a Curie- like 
1 /T behavior in both the spin-lattice relaxation rates and 
Knight shift at low temperatures. The STM experiments 
have revealed^ a sharp near-zero-bias peak around the Zn 
site, where the quasi-particle coherent peak, appearing 
near the superconductivity gap in the bulk case, is sup- 
pressed simultaneously. The zinc replacement has also 
shown a strong destructive effect on T c , and only sev- 
eral percentage of zinc doping can fully destroy the bulk 
superconductivity 7 . 

In order to explain the STM experiment beyond a con- 
ventional nonmagnetic impurity scattering treatment 8 , 
the existence of a local magnetic moment has been as- 
sumed, whose interaction with the d-wave nodal quasi- 
particles leads to a Kondo resonance peakSiiSiiiiiS. The- 
oretically it remains a great challenge to understand the 
microscopic origin of the magnetic moment found near a 
nonmagnetic Zn impurity. In one of recent attempts, it 
was interpreted^ as due to the binding of an S = 1/2 
nodal quasi-particle to the impurity, based on a mod- 
ified mean-field theory of the t — J model. The local 
staggered AF moments were also explained^ as a lo- 
cal spin-density-wave (SDW) ordering. But there still 
lacks a unified theory that can self-consistently explain 
the free moment and local AF ordering near the zinc im- 
purity without suffering a magnetic instability. Nonethe- 
less, various approaches have at least indicated that the 
zinc phenomena are quite different from those caused by 
nonmagnetic impurities in a conventional superconductor 
and have something directly to do with the nature of the 
underlying doped Mott insulators. A good understand- 
ing of the pure strongly correlated system, therefore, is 



quite essential in order to sensibly address the overall 
zinc impurity issue. 

In this paper, we approach the zinc problem by using 
a microscopic description of doped Mott insulators, in 
which AF correlations in spin degrees of freedom are sys- 
tematically described at various ranges as a function of 
doping concentration of holes. We show that an S = 1/2 
moment does emerge, naturally, near a zinc impurity, 
which is physically originated from an RVB pair in the 
original spin background. The latter becomes unpaired 
upon the Zn substitution, with one of its constituent spin 
being removed, together with the underlying Cu 2+ ion. 
In particular, we find that such an S = 1/2 moment 
cannot escape from the zinc impurity due a topological 
reason. It is a consequence of the nonlocal mutual en- 
tanglement between spin and charge degrees of freedom 
known as the phase string effect. Due to this effect, a 
neutral spin- 1/2 excitation (spinon) always carries a fic- 
titious 7T fluxoid as seen by the charge carriers and, what 
is more, even a vacancy on the lattice, such as a zinc 
impurity, also acts as a vortex when the charge carri- 
ers are condensed. Consequently, the condensate cannot 
survive in the presence of a zinc impurity unless a spin 
1/2, which carries an antivortex, is trapped nearby to 
compensate the vortex effect. 

In this bosonic RVB description, we show that the in- 
duced S = 1/2 moment around the zinc impurity dis- 
tributes like staggered AF moments due to the short- 
range AF correlations already present in the spin back- 
ground, which are frozen into a local AF ordering once 
the direction of the moment is fixed, say, by external 
magnetic fields. The calculated NMR spin-lattice relax- 
ation rates and uniform spin susceptibility are found to 
be in a systematic agreement with the experiments. The 
midgap spin excitations are also investigated, with the 
results consistent with the neutron experiment 1 ^. 

An important fact that we find is that all these novel 
properties are already exhibited in a zinc-doped state ob- 
tained by a 'sudden approximation': simply removing a 
spin sitting at the zinc site from the original pure system. 
A self-consistent adjustment of the RVB background be- 
yond such a sudden approximation only further strength- 
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ens the local trapping of the 5 = 1/2 moment as well as 
local AF staggered ordering around the impurity. There- 
fore, a Zn impurity in the high-T c cuprates constitutes a 
direct probe of the pure system according to the present 
theory. Namely, both the local 5 — 1/2 moment and the 
staggered AF moments 'induced' by a zinc substitution 
truthfully mirror the nature of the original doped Mott 
insulator. In this sense, the zinc substitution simply re- 
veals the 'secrets' of local physics already hidden in the 
zinc-free ground state. 

The remainder of the paper is organized as follows. In 
Sec. II, we present the topological reason that ensures an 
5=1/2 moment being trapped around a zinc impurity, 
based on the phase string effect in the bosonic RVB de- 
scription of the t — J model. Then in Sec. Ill, we present 
a generalized mean-field description for the doped Mott 
insulator with a zinc impurity. The detailed numerical 
results are given in Sec. IV. Finally, Sec. V is devoted 
to conclusions and discussion. 



II. TOPOLOGICAL ORIGIN OF THE LOCAL 
MOMENT AROUND A ZN IMPURITY 

A. Bosonic RVB framework 

We start with an effective description of the pure sys- 
tem of a doped Mott insulator, obtained 17,1 — based on 
the t — J model in an all-boson representation (i.e., the 
phase-string formalism) . 

The effective Hamiltonian is given by 17 ' 18 H e f / = Hh + 
H,, with 
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Here n\ a and nj 1 are spinon and holon number operators, 
respectively. By noting 0j(Z) = Im In (zj — z{) with Zj = 
%i + iy% representing the complex coordinate of a lattice 
site i, one has 
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and 4>1a is a uniform 7r-flux gauge field, satisfying 
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for each plaquette. It means that each holon behaves like 
a 7r fluxoid through A^-, which is felt by spinons in H s 
(here C denotes a closed loop and is the area encircled 
by it), and vice versa. Thus holons and spinons are mu- 
tually frustrated by each other nonlocally via the topo- 
logical fields, Afj and A^, which represents the phase 
string effect hidden in the t — J model. 

In this description, the superconducting state is real- 
ized by the holon condensation, (h') e lio / 0, while 
spinons remain RVB paired as characterized by A s ^ 0. 
Specifically, the superconducting order parameter can be 
written asi^iS 



A sc = A o^° 



+iz E I A iif + A E b ^ bia ~ N ^ ~ 6 ^ where A ° K ASh l and $S is defined by 



where h\ and bj denote the creation operators of bosonic 
holons and spinons, respectively, and Ay is the bosonic 
RVB order parameter determined either self-consistently 
or by minimizing the free energy. In the pure system, a 
uniform solution, Afj = A s , [ij S the nearest- neighbor 
(nn) sites], is usually obtained 18 . The Lagrangian multi- 
plier A in (J2J is introduced to enforce the global constraint 
of the total spinon number: 
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A unique feature in this model is the topological gauge 
field A'j, and A^, defined on a nn link (ij), as fol- 



which describes that each spinon carries a 2ir vortex in 
the phase of A sc (known as a spinon- vortesi*-) . Since 
spinons form singlet RVB pairs, these vortices and an- 
tivortices are generally cancelled out in (|ll|l such that 
the phase coherence of A sc can be established at low 
temperaturesi&iS. 

A non-superconducting state, with a finite pairing am- 
plitude A but is short of phase coherence, can be real- 
ized at higher temperatures where excited spinons dis- 
order the phase $f according to (jTT|) . Such a low- 
temperature pseudogap phase is called spontaneous vor- 
tex phase due to the presence of free spinon-vortices at 
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T c < T < T l 4*. The high-temperature pseudogap phase 
is defined at T v < T < T , where the holon condensation 
is gone such that the pairing amplitude A = 0, whereas 
the RVB order parameter A s still remains finite. The 
latter vanishes beyond To. 

The following discussion of local moments around zinc 
impurities will be mainly focused in the superconducting 
phase and spontaneous vortex phase, where the relation 
(Unf generally holds with A ^ 0. 



B. Topological origin of local moments 

Now let us consider a zinc impurity added to a pure 
system of the doped Mott insulator described by the 
bosonic RVB theory outlined above. 

In the high-T c cuprates, chemically it is a Cu 2+ ion 
that is replaced^ by a Zn 2+ . A Zn-potential may be con- 
sidered as a unitary potential 8 which pushes away both 
spin and charge from the zinc site. In the framework 
of the t — J model, a zinc impurity can be thus simply 
treated as an empty site with excluding the occupation 
of any electrons. 

Since a zinc impurity does not change the total charge 
of the system in the substitution of a Cu 2+ ion by a 
Zn 2+ , one may construct an effective theory by starting 
with that for a pure system and removing a neutral spin 
(spinon) from the system. A heuristic procedure is to 
imagine exciting a spinon at the would-be zinc site. With 
its spin being 'fixed', its exchange coupling with the sur- 
rounding spins is effectively cutoff. Nor a holon can hop 
to this site due to the no-double-occupancy constraint 
in the t — J model. Whether one removes or not such 
an isolated spin, the effective effect of a zinc impurity is 
created at such a site. 

Based on the bosonic RVB theory, corresponding to 
the creation of an isolated spinon at the would-be zinc 
site, a 27r vortex will then appear in the superconduct- 
ing order parameter (|1U|) via $ s or a tt vortex in the 
holon Hamiltonian via , in the superconducting 
or spontaneous vortex phase with the condensation of 
holons, (fv) ^ 0. In other words, each zinc impurity will 
always induce a nonlocal response (vortex current) from 
the charge condensate as shown in Fig. 1(a). 

Such a topological effect of a zinc impurity can be 
traced back to the nonlocal effect in the pure system 
of such a doped Mott insulator, known as the phase 
string effect: the motion of doped holes will always cre- 
ate string-like sign defects which cannot be 'repaired' at 
low energy in the t — J modeUi. Consequently spin and 
charge degrees of freedom are mutually 'entangled' in the 
pure system. An 'empty' (zinc) site then can be nonlo- 
cally perceived by the charge degrees of freedom to result 
in the above vortex-like response centered at this site. 

A more rigorous derivation is to start from the phase 
string representation^ of the t — J model with the pres- 
ence of an 'empty' site. Under the same RVB order pa- 
rameter A?., one obtains essentially the same effective 



Hamiltonians, QJ and J3J, except for that in (JSJ) the sum- 
mation does not include the 'empty' (zinc) site, denoted 
by io- Namely, 

4>i n = \ mo) - flj(io)] • 

Note that the summations in J2J and © should also ex- 
clude the zinc site, which are automatically ensured since 
both spinons and holons are not allowed at the site io, so 
that the definitions for Af, and A% remain the same as 
in the pure system. By noting 

for a closed loop around io, one finds that a zinc site is 
bound to an extra 7r-fluxoid which will always induce a 
vortex current in Q if the holons are condensed. Simi- 
larly, using the original definition of A in the effect 
of an empty site can be reexpressed by the replacement 
A° — > A°e I$i o , with <I> o giving rise to a 2ir vortex. Both 
are consistent with the previous argument based on freez- 
ing a spinon at the site i , which also lead to Fig. 1(a). 

Now it is natural to see why a zinc impurity will gener- 
ally induce a spin-1/2 around it in the superconducting 
and spontaneous vortex phases. In the superconduct- 
ing state, a Zn-vortex costs a logarithmically divergent 
energy and thus must be 'screened' by nucleating a neu- 
tral 5 = 1/2 spinon which carries an antivorte»i^ and 
is bound to the latter, as shown in Fig. 2b. Note that 
in the pure system, an isolated spinon excitation is not 
allowed in the superconducting bulk state due to the 
same reason, and only spinons in bound pairs (vortex- 
antivortex pairs) can be excited, which is known as the 
spinon confinements. It is also noted that in a differ- 
ent approach 20 , based on a similar all-boson formalism, 
an 'unscreened' current vortex is predicted around a zinc 
impurity if the S = 1/2 is trapped around which carries 
such a (anti)vortex. By contrast, such a (anti)vortex is 
always compensated in the present approach since a zinc 
impurity itself also induces a vortex as shown in Fig. 
1(a). 

In the spontaneous vortex phase, even though free vor- 
tices are thermally present in the bulk in a similar fashion 
as in Kosterlitz-Thouless transition, there still exists a 
logarithmic attraction between a Zn-vortex and a spinon- 
vortex at short-range, and a bound state (although not 
a confined state below T c ) between a Zn and an S = 1/2 
moment can be still present below T v . 

Therefore, there is a fundamental topological reason 
for an S = 1/2 moment to be trapped around a zinc 
impurity in the bosonic RVB theory of doped Mott in- 
sulators, in the superconducting and spontaneous vortex 
phases. 
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(b) 

FIG. 1: (a) A vacancy (zinc impurity) always induces 
a vortex-like supercurrent response in the superconducting 
phase due to the phase string effect, (b) To compensate such 
a vortex effect, a spinon, which carries an antivortex, has to 
be trapped around the zinc impurity, giving rise to a local 
S =1/2 moment. 



III. GENERALIZED MEAN-FIELD 
DESCRIPTION 

Once we have established the topological origin of the 
5 = 1/2 moment around a zinc impurity, a simple effec- 
tive description of the system with one zinc impurity can 
be developed. 

Note that after trapping a spinon nearby, the vortex 
induced by the zinc impurity is compensated by the an- 
tivortex carried by the spinon, as illustrated by Fig. 1(b). 
Then the system is topologically trivial at a distance suf- 
ficiently away from the impurity where the system re- 
mains the same as the pure system. The change of the 
state mainly occurs around the impurity with a charac- 
teristic scale comparable to the spin correlation length. 

We can construct a state based on the ground state 
I'J'o) of the pure system, by removing a local spin sitting 



at the site io, denoted by 

|* )zn = P io l*0> • (12) 

To leading order approximation, I'J'o) z n ma y ^ e regarded 
as a 'sudden approximation' of the true ground state 
l^ / )zn m t ne presence of a zinc. Both have the same spin 
and charge quantum numbers. By 'suddenly' removing a 
spinon at io, its original partner spinon in I'J'o) will be left 
around the site io in |vE , o)z n within the spin correlation 
length £. At distances larger than £, on the other hand, 
|\E'o)zn i s essentially the same as I'J'o) ■ Due to the topo- 
logical reason discussed in last section, the free spinon 
partner created in | v f , o)z n should remain trapped around 
the impurity. Therefore one expects |\E f o) Zn to have a 
good overlap with the true ground state |^)z n and can 
smoothly evolve into the latter under a weak local per- 
turbation. So it is reasonable for one to take (|12fl as a 
good variational form for the zinc problem. 

Since the holons are Bose-condensed in the supercon- 
ducting and spontaneous vortex phases which we are con- 
cerned, A^j in the spinon Hamiltonian @ can be simpli- 
fied as approximately describing a uniform flux with a 
strength 

per plaquette. 

Under the condition JT3J|, the spinon Hamiltonian H s 
in J2J can be straightforwardly diagonalized asi& 

H s = ^2 S m 7m Q 7ma + const., (14) 
by a Bogoliubov transformation 



with \u m \ = \J 2E "' an( l \ v m\ = y 2E m - ^ detailed 
treatment of this mean-field state in a self-consistent way 
can be found in Refsii2*2I. 

Based on the above mean-field description, the zinc- 
free RVB ground state I'J'o) is defined by "f ma I'J'o) = 0. 
Now we construct the states with one zinc being added to 
the system as discussed at the beginning of this section. 

Firstly, according to the sudden approximation, the 
trial state with a zinc at site io may be obtained by an- 
nihilating a bare spinon at site io with a spin index, say, 
—a: 

l*0) Zn = K-o |*0) ■ (16) 

Then, we allow the bosonic RVB order parameter Af,- to 
be adjustable around the zinc site to further minimize 
the total energy. Thus the final trial Zn state can be 
constructed in the following form 

|M>) Zn = Cb la _ a |*o[Ay) . (17) 
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Here C is a normalization constant and |\&o[Afj]) is the 
ground state of the Hamiltonian J5J under a fixed form 
of the RVB order parameter Af-, which will generally 
deviate from the uniform A s (of the pure system) around 
the zinc site. 

In order to address the dynamic and thermodynamical 
properties, one needs to further determine the elemen- 
tary excitations based on | v E') Zn defined in (|17fl . We take 
the following steps to make the construction. Firstly, by 
using the Bogoliubov transformation l|15|) , one has 

C6 io _ CT |*o) = C^iu ro<T (io)t> m 7,Lrl*o) 

m 

= &!*<>>, (18) 

with fL = C Y,m W ™,{io)v m r ll UJ and C = 

W2m \ w mcr (*o)| 2 "mj ■ Secondly, define a new class 

of single spinon creation operators j\ a as a linear com- 
bination of 7£ K7 's, which satisfies 

m m 

with FQ m = Cv m w mi j («o) such that /q ct is consis- 
tent with the definition in Q18)). A proper i* 1 can be 
then obtained by re-diagonalizing the Hamiltonian (I14|l : 
= Y, n jta,ot E nfX a fna + const., under a constraint 

foafoa = 1, with £?„ as the 'renormalized' spectrum 
in the presence of a zinc impurity. 

Then the ground state with a zinc is simply given by 
l*)zn = /ocr l^o) and the orthogonal (mean-field) exci- 
tation states are constructed by the creational operators 

fna ( n ^ 0) as follows 

IKa})zn= II (4a)" BB |*>Zn- (20) 

where v na denotes the occupation number at the state 
labelled by (n, a). 

Physically, in this approach is treated as a projec- 
tion operator, from the original ground state to the zinc- 
doped ground state, while f na (n ^ 0) creates spinon 
excitations which are ensured to be orthogonal to |4') Zn . 
The presence of a zinc impurity at site ig will be enforced 
by the constraint £ Q f^ a fo a = 1- 

Finally, it is noted that instead of treating Af ■ as un- 



Hamiltonian J5J will be kept at the zinc-free value to en- 
sure that the state remains the same at distances far away 
from the zinc site. The parameter po will be decided by 
enforcing the global constraint for the total spinon num- 
ber: 



^2(blb ia ) = N(l-S)-l 



(22) 



[compared to JSJl] in the presence of a zinc impurity with 
a given R.. In the following, most of results will be dis- 
cussed for the case of R — 1, i.e., within the sudden 
approximation, and then the stability of the results will 
be checked by allowing the variation of R. 



IV. PHYSICAL CONSEQUENCES 

Based on the above-constructed ground state and ex- 
cited states, we can straightforwardly calculate various 
related physical properties, in the presence of a zinc im- 
purity, by using the mean-field scheme similar to the pure 
system. The results are presented below. 

Local S = 1/2 moment. — As noted before, the ground 
state |\^) Zn defined in 117fl differs from the pure state 
|\&o) by a s P m 1/2- The change of the spinon distri- 
bution around the zinc in the ground state 1)17(1 is nu- 
merically determined, as illustrated in Fig. [3 where the 
hole concentration is fixed at 5 = 0.125, with the lat- 
tice size 16 x 16 and R = 1 chosen in Q21J1 . Compared 
to the zinc-free mean value 1 — 5 at a distance far away 
from the zinc site, the local density of spinons changes 
within a finite length scale near the zinc site as shown 
in Fig. |3 which accounts for the distribution of an un- 
paired spinon of S = 1/2 in 117|) created by the zinc 
substitution. Note that the spinon density is enhanced 
in the sublattice opposite to that of the zinc site, indi- 
cating that the unpaired spinon mainly stays there which 
reflects the fact that the underlying bosonic RVB pair- 
ing in l|17|l only involves spins at different sublatticcs 21 . 
It implies an AF spin configuration induced around the 
zinc as shown below. 

Staggered moments. — Corresponding to the above spa- 
tial distribution of the 5=1/2 moment, staggered (AF) 
moments are further shown in Fig. [3] around the zinc 
site. Note that the spin rotational symmetry breaking in 
Fig. [21 is because we fix —a =| in \^) Zn such that the 
total spin change upon a zinc substitution is AS Z = 1/2. 



restricted parameters in |\&o[Af,-]) 
simple site-dependence for Af • (ij E nn sites) in the fol 
lowing variational calculation, which is given by 



we shall assume a Tne state l*)z n in general is a spin doublet according to 



Af, = A s 

\ % J / nn 



l-(l-p )e 



(21) 



Here the bulk value A s is decided self-consistently in 
the zinc-free systemi2*2i. The zinc at site io will influ- 
ence (A?.-) within a radius R with »n determining the 

v t J / nn t ... 

strength. The Lagrangian multiplier A in the mean-field 



its definition in (JIT} (if |* ) 

is spin singlet). The length 
scale ^zn for the distribution of the local AF moments 
is essentially decided by the spin-spin corr elatio n length 
in the original RVB ground state, £ s « a y / 2/7r5( Refill). 
Fig. 0| shows the doping dependence of the scale ^zn for 
the distribution of the local spin moments, defined by fit- 
ting (Sf) ~ (-l) l 5 exp(-|z - i | 2 /£l„). Indeed we find 
£ Zn — £«• For a hole doping S = 0.125, £ Z n — 2.22a. For 
a typical zinc doping in experiment, Sz n = 0.03, the aver- 
age distance between two zincs is dz n ~ <V ' \/5z n — 7.07a, 
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FIG. 2: Spinon density distribution (n\) = \Y,^ b L b ' 
around the zinc impurity. The lattice is 16 x 16 with dop 
ing 5 = 0.125. 
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FIG. 3: The distribution of (Sf) near the zinc impurity, with 
the scan along the dashed direction shown in the inset, where 
the zinc site is marked by the filled circle. 



which is larger than £zn at 6 — 0.125. In this dilute case, 
the correlation effect among different zinc impurities can 
be neglected and one may only focus on the single zinc 
effect. 

Low-energy spin excitations. — One can further exam- 
ine the dynamic properties of the induced local AF mo- 
ments. One physical quantity is the spin-lattice relax- 
ation rate of 63 Cu measured in NMR/NQR experiments, 
which is decided by the imaginary part of the spin corre- 
lation function % (g, cj) as follows 



63 TlT 



E 



A 2 (q) Imx(q,u N ) 



LON 



(23) 



Here the structure factor A (q) = A + 
2B (cos q x a + cos q y a) and A — —AB when the ap- 
plied magnetic field is orthogonal to the Cu02 plane 1 . 
So the main contribution will come from the AF 
correlations near the AF momentum (tt, 7r). Fig. [S] 
shows the theoretical calculations. In the pure system 
of the bosonic RVB state, a pseudogap opens up in 
the spin excitations, resulting a suppression of l/ 63 TiT 



FIG. 4: The size of the spatial distribution of the local 
moments, ^„ (solid), and the bulk spin correlation length, 
£ s = ^/2/<$7r (dotted), vs. doping 8. 



at low temperature (solid line with crosses in Fig. EJ). 
However, for the zinc-doped system, the 'spin gap' 
shown in Fig. [SJis filled up by a Curie- type contribution, 
l/ 63 TiT oc l/T, due to the presence of a free moment 
with a staggered distribution around the zinc site. The 
spatial distribution of l/ 63 TiT is given in Fig. [f)|at a low 
temperature (T — 0.0067 J), which clearly shows that 
the Curie-type signals are located near the zinc site, 
with its maximum at the nn sites of the zinc, consistent 
with the distributions of the staggered moments in Fig. 

El 

It is noted that the calculation is performed at a 16 x 16 
lattice with only one zinc, with an effective zinc doping 
concentration equal to 1/256 ~ 0.004. In the realistic 
case with more zinc impurities, the intensity of l/ 63 TiT 
on average is expected to be proportional to the density of 
zinc impurities at low temperatures, so long as the zinc 
concentration is not too large such that the single-zinc 
approximation used above is still valid. Thus one should 
multiply a factor 5 ~ 8 to the low-temperature part of 
l/ 63 TiT for the averaged case in Fig. [S] (solid curve with 
squares) in order to compare with the experimental case 
with <5 Z „ = 0.02 - 0.03. 

Furthermore, the uniform spin susceptibility can be di- 
rectly calculated by Xu oc Rex (<7 = 0, u — > 0).. \u also 
show a pseudogap behavior at low temperatures in the 
pure state, which is replaced by the Curie- like l/T be- 
havior near the zinc site due to the contribution from the 
local moment as shown in Fig. \7\ 

Dynamic spin susceptibility. — The imaginary part of 
the dynamic spin susceptibility, Imx(g,w), can be di- 
rectly measured by inelastic neutron scattering. In the 
bosonic RVB mean field theory, a resonance peak at the 
AF wave vector Qaf = (tTj tt) is present in the dynamic 
spin correlation function in the superconducting phase. 
For the hole doping 5 = 0.125, we obtain the resonance 
energy E g » 0.53J2I. Upon the zinc doping, the energy 
Eg of the resonance peak has changed little (see Fig. [SJ) 
in the bulk. But a zinc impurity does induce some new 
states at lower energies as shown in Fig. [SJ which reflects 
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Temperature (J) frequency m (J) 



FIG. 5: Contributions to l/ 63 TiT from different sites are cal- 
culated. Solid curve with crosses: from the site far from the 
zinc impurity; Dashed curve with triangles: the nn site near 
the zinc; Dashed curve with circles: the next nearest neigh- 
bor (nnn) site near the zinc; Solid curve with squares: average 
over all sites in a 16 x 16 lattice with one zinc (see text). 



FIG. 8: Dynamic spin susceptibility at AF wavevector Qaf = 
(7r,7r) with 8 = 0.125. The high-energy resonancelike peak 
at Eg ~ 0.53J is from the bulk (zinc free) state, while the 
low-energy excitations indicated by the arrow are new ones 
induced by the zinc impurity. 
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FIG. 6: Distribution of the contributions to l/ 63 TiT from 
individual sites near the zinc impurity, at temperature T = 
0.0067J. 
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FIG. 7: Uniform spin susceptibility in the pure system is 
shown by the solid curve with full circles; The case with one 
zinc is illustrated by the solid curve with triangles; The dashed 
curve is a fit by x = 0.2390/T. 



the modified spin excitation spectrum near the zinc, ac- 
companying the emergence of a local moment. Note that 
the weight of such zinc-induced modes in Fig. [S] should 
be enhanced with a finite concentration of the zincs. 

Effect from the holon redistribution. — So far the mean- 
field results are obtained based on the assumption that 
the Bose condensed holons are uniformly distributed in 
space. But in the presence of a zinc, the holon density 
near the zinc impurity should be generally suppressed, 
to be consistent with that the spinon density increases 
around the zinc site due to the no double occupancy con- 
straint. In the bosonic RVB theory, since the holons will 
influence the spinon part by the gauge field A'^ defined 
in (jS)) , the suppression of the holon density around a zinc 
will cause additional effects on the local spin dynamics, 
which is considered below. 

As shown in Fig. [51 the intensity of the zinc-induced 
low-energy spin excitations seen in Fig. [S] will be en- 
hanced with its energy scale further reduced if a holon 
density reduction is taken into account with a profile 
given by the inset. Meantime, the local staggered mo- 
ments will also be increased (Fig. HOfl under the same 
holon distribution. Although the above calculations are 
not based on a self-consistent scheme, which is generally 
quite difficult, all the important features found previ- 
ously in the uniform profile of the holon distribution are 
kept unchanged qualitatively, except for that the anoma- 
lies are further strengthened due to the reduction of the 
holons which makes the area near the zinc closer to the 
half-filling. 

Local stability of the mean-field theory. — The stability 
of the results discussed above based on the conjectured 
ground state l|17|) can be further examined by tuning the 
variational order parameter defined in l|21|l . By changing 
the parameter R, the range of the zinc effect on the RVB 
parameter (Af^) in l|21l) can be continuously adjusted. 

As shown in Fig. 1111 the superexchange energy (-ffj)zn 
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FIG. 9: Dynamic spin susceptibility at Qaf = (tt, 71") is shown 
as the solid curve when the holon density is suppressed locally 
around the zinc (see the inset). For comparison, the dashed 
curve illustrates the case for a uniform holon density distri- 
bution. 




lattice site 



FIG. 10: The distribution of |(5T)|, corresponding to the pro- 
file of the holon density shown in the inset of Fig. [5] is plotted 
as the solid curve, as compared to the dashed curve for the 
uniform holon distribution. 

is found to monotonically grow with the increase of R for 
both uniform and non-uniform profiles of the holon distri- 
bution, indicating that the sudden approximation state 
|^o)zn — \^(R = l))zn remains locally stable, which pre- 
vents the local 5 = 1/2 moment from leaking far away 
from the zinc impurity When the holon redistribution 
discussed above is considered (see the dashed curve in 
Fig. Ill|) , this local stability is further strengthened, since 
the frustration effect on the spin dynamics, which comes 
from the holon motion, is weakened around the zinc im- 
purity due to the reduction of the holon density. 



V. CONCLUSION 

In this paper, we have developed a microscopic de- 
scription of the zinc doping effect in the cuprate super- 
conductors based on an effective theory, i.e., the bosonic 



FIG. 11: The average superexchange energy (Hj)/(2Nl) per 
bond (Nl denotes the total number of lattice sites excluding 
the zinc sites) as a function of the parameter R in \Y2 \t . Solid 
curve: uniform holon distribution. Dashed curve: locally sup- 
pressed holon density around the zinc indicated by the inset 
of Fig. 10. 



RVB theory, of the t — J model. The unusual effects of 
the zinc doping observed in the experiments have been 
explained as a direct consequence of the mutual nonlo- 
cal entanglement between the spin and charge degrees of 
freedom (i.e., the phase string effect) in the doped Mott 
insulator. Once a local moment is formed due to such 
a topological mechanism, the effect of zinc doping can 
be simply described by a 'sudden approximation', which 
directly translates the short-range RVB pairing, present 
in the spin background state of the pure system, to the 
AF spin distribution around a zinc impurity, once the 
latter is introduced to the system. The NMR spin relax- 
ation rates, uniform spin susceptibility, and the induced 
low-lying spin excitations, obtained within the sudden 
approximation and beyond, have consistently painted a 
unified picture of spin correlations near the zinc, which 
is shown to be stable locally and in an overall agreement 
with the experimental observations. In particular, the 
theory predicts that the range of the distribution for the 
local moment is inversely proportional to the square root 
of doping concentration. 

In this work, an important property in the zinc prob- 
lem has not been discussed so far. That is the behavior 
of quasi-particle excitations and the single-electron tun- 
nelling properties. In the bosonic RVB theory, an elec- 
tron is composed of a holon and a spinon, together with 
a nonlocal phase string factor—. In the superconduct- 
ing phase, it has been shownS^ that a quasi-particle is 
stable due to a confinement of these holon-spinon and 
phase-string objects. So the overall low-energy single- 
particle spectral function is expected to be the same as 
in a d-wave BCS theory. However, the composite struc- 
ture is predicted to be seen at high-energies, which may 
explain the 'coherent peak' in the antinodal regime in 
the pure system as due to the spinon excitation— In the 
zinc doping case, due to the trapping of a free moment 
(spinon) near the zinc, the high-energy spinon contribu- 
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tion in the single-particle spectral function can become a 
zero-biased mode. Detailed study of this property needs 
the knowledge about the charge degree of freedom as well 
as the confinement in the superconducting phase, which 
is beyond the scope of this work. This issue and other 
problems, like how the superconducting transition tem- 
perature will be affected by the local moments, will be 
further investigated within the bosonic RVB description 
in future. 
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